## This code creates Figure 1.

## This code requires Google Mpas API key to run. Please run register_google([Insert Google API Key]).

# Figure 1
 
mm_basemap <- get_googlemap("Yangon, Myanmar", zoom = 5, maptype = "terrain")

bago <- spTransform(region,crs(konbaung))

ggmap(mm_basemap) +
geom_polygon(aes(x=long, y=lat, group=group),color='black',fill='white', data=konbaung,linetype=3,size=0.5,alpha=0.8) +
geom_polygon(aes(x=long, y=lat, group=group),color='black',fill=NA, data=myanmar,size=0.5) +
geom_polygon(aes(x=long, y=lat),color='black',fill='black', data=bago,size=0.5)



